function xdeltaNSslopeFn = deltaNSslopeFn(x)
xdeltaNSslopeFn = (-1250*((1041*x(1))/1250 + (1041*x(2))/1250 + (1041*x(3))/1250 -  ...
   (1903*x(4))/5000))/1041 ...
;